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Abstract 

We have shown that it is possible to formulate the consistent and probability-preserving de¬ 
scription of the CP-symmetry-violating evolution of a system of decaying particles. This has been 
done within the framework of quantum mechanics of open systems. This approach allows the 
description of both the exponential decay and flavour oscillations. We have solved explicitly the 
Kossakowski-Lindblad master equation for a system of particles with violated CP symmetry and 
found the evolution of any observable bilinear in creation and annihilation operators. The choice 
of a concrete observable can be done by the proper choice of initial conditions for the system of 
differential equations. We have calculated the evolution as well as mean values of the observables 
most interesting from the physical point of view, and we have found their lowest order difference 
with the CP-preserved values. 
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I. INTRODUCTION 


In recent years Bell ineqnalities [l| were tested in systems of correlated nentral K 0,3 
or B mesons 0. However, the quantum mechanical analysis of such systems encounters 
an important difficulty, namely, the irreversibility of time-evolution of unstable particles. 
Only the complete system consisting of the decaying particle as well as the decay products 
undergoes unitary evolution, which is actually described by quantum held theory. But in 
Einstein-Podolsky-Rosen (EPR) correlation experiments 0, it is more useful to consider the 
decaying particles only, neglecting the evolution of decay products. In the usual approach one 
introduces a non-Hermitian Hamiltonian, as it was done in the classical works of Weisskopf 
and Wigner 0. However, the non-Hermitian Hamiltonian does not provid an unambiguous 
way of calculating the probability of Ending the system consisting of a few such particles in a 
fven state after the measurement (such an approach causes also other ambiguities, see e.g. 

for a discussion). Indeed, using this approach, we must reinterpret quantum mechanical 
results with the help of probability theory, especially if we would like to analyze correlation 
experiments with unstable particles, namely, neutral K or B mesons 0 . From the quantum 
mechanical point of view, this means that unit trace and positivity of the density matrix for 
the system under consideration must be preserved 0. Fortunately, it is ^ssible to resolve 
the above-mentioned problems using the open quantum systems theory The idea that 
unstable particles can be treated as open quantum systems was proposed Erst by R. Alicki 
m and was developed by various authors in difierent contexts (see e.g. [H, 13|). Recently, 
this approach was successfully applied to the systems of particles with Eavour oscillations 
(like in the case of neutral kaons or B mesons) [14-10, and it has been used successfully in 


the description of EPR correlations and evolution of entanglement in K°K° system 17|, [18 


In this theory the evolution of a quantum system is described by the master equation [19 


which can be treated as the replacement for either von Neumann or Heisenberg equations, 
for the picture of quantum mechanics in use. Here, we follow the approach presented in 
T^. where systems with an arbitrary number of particles were described with the use of the 


second quantization formalism, which seems to be the most natural language for a system 
with varying number of particles. 

The paper is organized as follows. In Sect. [ITlwe introduce the notation and conventions 
used through the paper. In Sect. [TTIl we introduce the master equation for a system of 
neutral kaons with violated CP symmetry and we End the evolution of the observables in 
the Heisenberg picture, and, in Sect. HVl we analyze the evolution of the total number of 
particles and strangeness observable as well as the numbers of each of neutral kaon Eavor. 
Finally, we conclude the paper in Sec. |V| 


II. PRELIMINARIA 


We start with introducing some notations, conventions and deEnitions used throughout 
the paper, as well as give experimental values of some parameters important in the descrip¬ 
tion of the systems of K° and K°. 

We deEne one-particle states |K°) and |K°) as 

|K°) = at |0) , 

|K°) = fet |0) , 
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where a and b fulfills the usual canonical commutation relations 


[a, a^] = [b, 6^] = 1, (2a) 

[a, b] = [a^, 6^] = 0 . (2b) 

States |K°) and |K°) are orthonormal, i.e. 

{K°|K°) = (K°|K°) = 1, (3a) 

(K°|K°) = 0 , (3b) 

and they are eigenstates of the strangeness operator S = a^a — Wb\ 

^ |K°) = |K°) , (4a) 

S |K0) = - |K0) . (4b) 


For our convenience we will use the following short-hand notation for multiparticle states: 

|#K° = n,#K° = h) = |n,n) . (5) 

Therefore 

| 0 ) ( 6 ) 

Vn!n! 

and 

S |n, n) = {n — n) |n, fi) . (7) 

However, the states |K°) and |K°) are not eigenstates of time evolution. Actually, under 
time evolution the “well behaving” states are |Kg) and |Kl), defined as follows 


|K°)=p|K°) + g|K0) , (8a) 

K°)=p|K°)-g|K°) , (8b) 


with mean life times ts = 0.8954 x 10 ^‘^s and T£ = 5.116 x 10 ^s, respectively j^ . 
Moreover, li p ^ q then CP symmetry is violated. Indeed, 


p _ 1 + e 
q 1 — e 


(9) 


with 


ipr+kr = i, (loa) 

ipr-kr=^L, (lob) 


and Al = 2 Re(e)/(I + |ep) is a measure of violation of CP symmetry (hereafter we use the 
so-called Wu-Yang phase convention 


are Al = 0.332% and 
since 


el = 2.228x 10 


21 

in 

-3 

20 


which e = e). Experimentally obtained values 
Notice also that the basis (jS]) is not orthogonal 


(K°|K°) = A^. 


( 11 ) 
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III. MASTER EQUATION 


The aim of this section is to hnd the evolution of the multiparticle system of and 
K°. We know that the evolution of a single neutral kaon can be described as an open 
system obeying the Kossakowski-Lindblad master equation [l^. In this case the one-particle 
Hamiltonian for this master equation was written in the base ([8]) and reads 


H = {™S IKS (KS + rUL |K») (K"! 

-A, [{m - jAF) |K») (K»| + (m + JAF) |K«) {K»|] } , (12a) 

while Lindbladians are of the form 


Li = 

L2 = 


1 

1 

l-Al 


(|0) (K»| - A, |0) {K»|) , 

|0>(K»|-/i,(vrr 


r — iAm'^ 

WrUj 


|0) (Kg 


(12b) 

(12c) 


and K = —{L\Li + L\L2)/2, where m = (rn^-|-mL)/2 = 497.614 MeV/c^ is the K° mean 
mass, Am = mi — mg = 0.5293 x 10^° h/s [^, T = (Tg -|- rL)/2, and AT = T^ — Ty,. 

The state of a single particle system under the evolution given by the master equation 
with operators flT^ is in general a mixed state of a single particle and vacuum (see fl^b 
After the projection of this state on the single-particle sector one gets the state obtained by 
the non-trace-preserving equations, considered usually in the literature (see e.g. jsl). 


A. Master equation in second quantization formalism 


Now, let us go to the description of a multiparticle system. We achieve this by the 
replacement of intertwining operators with the appropriate combinations of annihilation 
and creation operators, as was done in [l^. In the base ([T]) this Hamiltonian and these 
Lindbladians take the form 


H = m (a'a 


(a^a -|- b^b) — ^ (Am -|- ^AgAV^ a^b — ^ (Am — lA^AT) (13a) 

1 - /\ T 1 - J\. T 


qp 


and 


L, = Jrs-Al 


-|- Am^ / p 


1 + 4 


-Q 


1- A; 


b] , 


L, = 


P 


l + A, 


v^ + 4 


r — iAm 




vri J i-2i 


. T — iAm\ , 
Al -4 


c/Fi y 


(13b) 

(13c) 


K = —It (a4 -|- b^b'j - ^ ~ iAgAm) a^b - ^ (1'^^ + i^iAm') b'^a 

(13d) 


qp 


Let us consider a single particle state described by a density matrix of the form 
p = Py |K°) (K°| + p 2 |K°) (K°| + w |K°) (K°| + w* |K°) (K°| + (1 - pi - ps) |0) (0| (14) 
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with 0 < Pi < 1, 0 < p 2 < 1, 0 < Pi + p 2 < 1, and < piP 2 - The evolution of this state 
in the Schrodinger picture is governed by the master equation 

dtp{t) = -i[H, p{t)] + {K, p{t)} + Lip{t)L\. (15) 


By a straightforward calculation, using relation (|H]) between the bases {|Kg), |Kl)} and 
{|K°), |K°)}, one can check that the choices of Eqs. flT^ and flT^ for operators appearing 
in Eq. f[T5]l lead to the same evolution equation for p{t). 

Moreover, the evolution of the two particle system obtained from Eqs. 0131) and 015p lead 
to the evolution of the two particle state obtained by taking a symmetrized (the particles 
are indistinguishable) tensor product of single particle evolution (see 0). 

Now, consider the Kossakowski-Lindblad master equation for the evolution of observable 
VL{t) in the Heisenberg picture 

dtn{t) = C[Q{t )], (16) 

where 

+ + . (17) 

i 

and we choose H and Li {i = 1, 2) in the form ([13]). We assume that Q{t) can be written as 
a bilinear form in annihilation and creation operators 

^{t) = uJaa{t)a^a + uJab{t)a'^b + uJba{t)b'^a + ujbb{t)b^b , (18) 

with the condition Ubaif) = guarantee that r2(f) is Hermitian. 


B. General solution 

Here we want to hnd the evolution of any observable fl{t) which is bilinear in annihilation 
and creation operators. To achieve this aim, we begin with the following observation; 

C[a^a] = —Ta^a — ^ — iAm) a)b — ^ (1"^^ + iAm) b^a , 

C[a^b] = —Ta^b — ^ ^^ — iAm) a^a — ^ ^ + iAm) b^b , 

C[b^a] = -Tb^a -(i AT + iAm) a^a - {^AT - iAm) b^b , 

£[6^6] = —Tb'^b — ^ ^ (1^^ + iAm) a^b — ^ — iAm) b'^a . 

We see that if Vt{t) is of the form 0181) at a given moment of time Iq, it must preserve this 
form for all the time f > to¬ 
ff we take into account the linearity of Q{t) and linear independence of operators a^a, 


(19a) 

(19b) 

(19c) 

(19d) 
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a^b, b^a, and 6^6, we get the following system of first-order differential eqnations for lo: 


uJaait) = -TuJaait) - (|Ar - iAm) LJabit) - (^Ah iAm) LJbait) , (20a) 

t + t + 

(^ab{t) = - (|Ar - iAm) Uaa{t) - Tuab{t) - (iAh iAm) Ubb{t) , (20b) 

J- — 

uJbait) = - (|Ar -h iAm) Uaait) - Tubait) - (^Ah - iAm) Ubbit) , (20c) 

i — Al i — Al 

Abb{t) = - (|Ar iAm) ujab{t) - (^Ah - iAm) Uba{t) - Tujbb{t). (20d) 

J- — -L — 

Using straightforward methods we can hnd that the general solntion of (120]) is of the form 

g-rt 

UJaa{t) = -Y- { [cosh (^AH) -h COs(Amf)] Uaaifi) 

— - [sinh (iAFf) — isin(Amf)] a;ab(0) — --[sinh (iAFf) -|- isin(Amf)] a;ba(0) 

P 1 + Al q 

1 — At 1 

-— [cosh (iAFt) — cos(Amf)l (^^^(O) > , (21a) 

1 + J 

—Ft ( 

(^ab{t) = ^-Y [®i^^ (i^rt) - isin(Amf)] ^^^(O) 

-f [cosh (|AFt) cos(Amt)] Ci;a6(0) ^ [cosh (|AFf) — cos(Amt)l Uba{0) 

I + Al \qJ 

- -[sinh (^AFt) isin(Amf)l ^^^(O) 1 , (21b) 

^ + AlQ I 


(^ba{t) 


e 


-rt 


2 


+ Al ( 

1 - V 


- [sinh (iAFf) -|- isin(Amf)l UaaiO) 

1 — AlP 

[cosh (|AFt) — cos(Amt)] Ci;afe(0) + [cosh (^AFt) -|- cos(Amt)] a;fea(0) 

— - [sinh (iAFf) — i sin(Amt)] a;fcfe(0) j- , (21c) 


(^bb{t) = 


,-rt 


1 + g 

1- AlP 


I “ cos(Amf)] a;„a(0) 

I 1 — 

[sinh (iAFf) -|- isin(Amf)] a;afe(0) — - [sinh (|AFf) 

-|- [cosh (|AFf) -|- cos 


— isin(Amf)] a;ba(0) 
(Amt)] (nbb(O)} . (21d) 
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The Hermicity of fl{t) is preserved (i.e. the condition Ci;ab(0) = Ci;^a(0) implies ojabif) = 
for all the time f > 0), since from Eqs. ([2]) and fllUp we have 


/pV _ l + Ai q 

\q) 1 -Alp' 


( 22 ) 


Now let us observe that the different choices of initial conditions for ca’s allow us to get 
the time evolution of different physically interesting observables of the form flTSD . 


IV. EVOLUTION OF THE OBSERVABLES AND THEIR AVERAGES 

Now we are prepared to hnd the evolution of some physically interesting observables, 
as well as the time dependence of their mean values. We begin with the total number of 
particles and strangeness observables, next we analyze the numbers of each neutral kaon 
flavor. 


A. Total number of particles and strangeness 

For the total number of particles observable N, we have iV(0) = a^a + b^b, so u)aa{0) = 
= 1 and a;afe(0) = a;fea(0) = 0, and hnally we get 

N{t) = I ^ ^ ^ [cosh (|Art) + Aicos{Amt)~\ a^a 

— ^ ^ ^ - [sinh (IArt) — zAisin(Amt)] a^b 

1 Q 

--——— [sinh (|Art) + zAi sin(Amt)] b^a 

+ - -— [cosh (IArt) — Aicos^Amt)] b^b\ . (23) 

I — Al J 

The expectation value of iV(t) in the state \n,n) is 

g-rt 

(A(t)) = { [cosh (|Art) — A\ cos(Amt)] {n + n) 

—Al [cosh (IArt) — cos(Amt)] (n — h)} (24) 

and is depicted in Fig. [T] If we take = 0 in Fq. fl24D . we obtain the case with preserved 
CP symmetry, for which we recover the usual result for mean number of particles 

(^(^))cp = I (n + n). (25) 

The leading part of the difference between UP-violated and UP-preserved values is 

(A(t)) - {N{t))(^p = -Al [| cos(Amt)] (n - n) + O(A^), (26) 

and is shown in Fig. (H too. 
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FIG. 1 (Color online). Mean number of particles vs time for n + n = 1,... ,5 (bottom-up) and 
the difference between {N{t)) for CP-violated and CP-preserved case vs time for n -|- h = 4 and 
n — h = 0,...,4 (top-down). 


For strangeness S we have S'(O) = a'^a — b^b, so a;(ja(0) = —a;6fe(0) = 1, ^^^(O) = ^^^(O) = 0, 


and 
S{t) = 


l + Ai 


\Al cosh (IArt) -f cos(Amt)] a^a 


1 p 


1 + Aiq 

1 q 

I- Alp 


[A^sinh (iAFt) — isin(Amt)] a}b 
\Al sinh (|Art) -b i sin(Amt)] b^a 


1-Al 

The expectation valne of S{t) in the state |n, h) is therefore 


[A^^cosh (iAFt) — cos(Amt)] 6^6 > . (27) 


{S{t)) = -—^ I [cos(Amt) — A\ cosh (^AFt)] (n — n) 

+Al [cosh (^AFt) — cos(Amt)] {n + n)} (28) 


and is shown in Fig. [2J As previonsly, when A^ = 0 we get the CP-preserved valne 





















f(ns) 

FIG. 2 (Color online). Mean strangeness vs time for n + n = 4 and n — n = 0,..., 4 (bottom-up) 
with the strangeness oscillation phenomenon explicitly visible when n — n ^ 0 and the difference 
between {S{t)) for GP-violated and CP-preserved cases vs time for n -|- h = 1,..., 5 and n — n = 0 
(bottom-up). 


{S{t))(jp = e cos{Amt){n — n), 


(29) 


SO the leading term of the difference between CP-violated and CP-preserved valnes is 


{S{t)) - (S'(t))cp = Al (e -b e - e cos(Amt)] (n + n) + 0{A‘l) (30) 


and is shown also in Fig. |5J 
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B. Number of K° and K 


"0 


For the numbers of and K° we have A^ko(O) = and iVj^o(O) = 6^6, so 


-,-rt 


iVKo(t) = —— < [cosh (|Arf) + cos(Amf)] a^a 


P 


[sinh (|Arf) — isin(Amt)] a^b 


^ [sinh (iAFf) + isin(Amf)] b^a 


I- Alp 


+ '\—^ [cosh (iAFf) — cos(Amt)] b^b 
1 — A.r, 


(31a) 


and 


A'ko(«) = 


,-rt 


1-dl, 

1 + Ai 


[cosh (IArt) — cos(Amt)] a^a 

[sinh (|Art) + isin(Amt)] a^b 


1 + AlQ 

q 


— [sinh (|Art) — tsin(Amt)] 6^a 


[cosh (^AFt) + cos(Amt)] b^b > . (31b) 


Their mean values are 


g-rt r _l_ 

(AKo(t)) = < [cosh (^AFt) + cos(Amt)] n + -—— [cosh (|AFt) — cos(Amt)] n 

(32a 


and 


-,-rt 


(%o(t)) = 


1-dlz 

1 + Ar 


[cosh (iAFt) — cos(Amt)] n + [cosh (^AFt) + cos(Amt)] n | , 

(32b) 

respectively. We show these time evolution in the Figs. OandlH Notice that after a suitable 
period of time (approx. 8 ns) different initial states can give the same mean number of 
or K°. 

If we put Ai = 0, we get again (FP-preserved values 


and 


(AKo(t))cP = 


(^ko(^))cp - 


j-Tgt _|_ p-Pit , 

y 2 

,-rst _j_ g-r,Lt 77,-1-72 


+ e “cos(Amt)- 


e cos(Amt) 


n — n 


2 2 ' ' 2 
The leading term of the differences for the CP-violated and CP-preserved vanes are 


(33a) 

(33b) 


(AKo(t)) - (lVKo(t))cp = Al [| (e -1- e - e cos(Amt)] n + 0{A\) , (34a) 

(AKo(t)) - (A'ko(^))cp = -Al [^ cos(Amt)] n -f O(A^). (34b) 
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FIG. 3 (Color online). Mean number of K® vs time for n = 1,2,3 (bundles bottom-up) and 
n = 0,1, 2 (bottom-up inside bundles). 


C. Number of Kg and 

Let us return to the short and long living states of neutral kaon, |Kg) and |Kl), re¬ 
spectively. If we define another two pairs of annihilation and creation operators, such that 


|Kg) = 4 |0) . (35a) 

IKS’) = 4 |0) , (35b) 

we can hnd, with use of Eq. ([S]), that 

cs = p*a + q*h, (36a) 

Cl = p*a — q*h. (36b) 

The operators cs and cl fulhll almost the usual canonical commutation relations with the 
following exception: 

[cs, 4] = [cl, 4] = ^L • (37) 

This reflects the fact that |Kg) and |Kl) are not orthogonal, see Eq. ffTT]) . 
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FIG. 4 (Color online). Mean number of vs time for n = 1, 2,3 (bottom-up inside bundles) and 
n = 0,1, 2 (bundles bottom-up). 


Now we define multiparticle states for Kg and Kl as 


(J '\n 

|#K° = n)^|ns) = ^|0) , 

Vnl 

(J )n 

mi = n)^\n,)=^-jL\0) . 

Vnl 


(38a) 

(38b) 


These states can be expressed in terms of the sates \n, k) as follows 


Wl) = I” “ 


(39a) 

(39b) 
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For the number of Kg, iVK;o(0) = c\cs = + ^—^b^b, so 

1 


,-ri 


NKoAt) = 


2 I 1 + 

1 p 


e 2 ^^* + + 2Al cos(Amt) 


a^a 


+ 


1 + g 

1 q 


1- Alp 


e 2 ^^^* - A\e 2 ^^^ + sin(Amt) 


e 2 ^^* - A|e 2 ^^* - sin(Amt) 


a^6 
6'^ a 


1-A, 


e 2 ^^"* + A^e 2 ^^* - 2Ai cos(Amt) 


6^6 I . (40) 


Similarly, for the number ofK^, Ako(O) = ^ ^'^ 2 ^'^ pb^(^+ ^ 2 ^ b'^b, 

and 


,-rt 


AKo(t) = 


2 11 + Ai 

1 p 

1 + AlQ 
1 q 


+ A^e 2 ^’"* + 2Ai cos(Amt) 


a'a 


1- AlP 


e 2 ^’"* - A^e 2 ^^^* - 2iAL sin(Amt) 


- A^e 2 ^^* + 2iAi sin (Amt) 


a^b 

b^a 


1-A, 


+ A^e 2 ^^^* - 2Az, cos(Amt) 


b^b > . (41) 


Now we are prepared to hnd the expectation values of iVKo(t) and iVKo(t) in the states |n 5 ) 
and Iul), respectively. 

First, observe that 


a)b n 5 ) = -a)a n^) , 

p 

a^b Iul) = 

--a^a Ini) , 

p 

(42a) 

b^a \ns) = -b^b n^) , 

b^a Iul) = 

-h^blriL) . 

(42b) 


Thus 


A^Ko(f) \ns) = 


-.—Ft 




e-iAri _ ^2 j<An>A ;v(0) - Al (e"2^''' - e‘^“‘) S(0)] \ns) (43a) 


and 




a-rt 


1-^i 


e^Art _ ^2^^-iAmt 


A(0)-Ai e2 


^^AFt _ ^—iAmt 


Of course |ns) and Iul) are eigenvectors of A(0): 


^( 0 )J Ini) . 

(43b) 


A^(0) \ns) =n\ns) , 
A^(0) Iul) = n \nL) , 


(44a) 

(44b) 
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but they are not eigenvectors of S'(O) since 


>^(0) = n \ns) - 2 ^ V \n -k,k) , 

-5(0) In^) =n\nL) - 2'^{-l)’^ki ( 

k=0 V 


pn-kqk\n_k^k) . 


(45a) 

(45b) 


Notice that in view of Eq. fl45l) the states |n 5 ) and Iul) are not eigenvectors of N^o and 
A^ko respectively, as one can naively expect. However, the expectation values of *5(0) in the 
states |?7,5) and {ul) are well established and can be calculated as 


{ns\S{0)\ns) = Alu, (46a) 

(nil 5(0) Irii) = HiU, (46b) 


because 

E*-(!)ipi^‘"-‘'i9i“ = "i9t (ipp + = "kr (V) 

k=0 ^ ' 

In conclusion, we get that the mean number of Kg and K® in the states In^) and In^), 
respectively, evolves in time according to the Geiger-Nutall law 


(n^l N^{t) \ns) = ne , (48a) 

(nil N'ko (^) I’^l) = . (48b) 


The states jn^) and jn^) are not orthogonal, thus we expect nonzero expectation values 
of N^o(t) and iV^o in the states jn^) and jn^), respectively. Indeed, 


'-L 

e~^^Ai 


^Kg(t) \nL) = L 


—iAmt 


S{0) - Al [e 


—iAmt Art 


62 ^“ Ar(0) |ni) (49a) 


and 


^K?(^) \^s) = 


e ^^Ai 

L 


JAmt _ j^2 g—iArt 


5(0)-Hi e 


AAmt 


— e 2 


-iArt 


thus. 


{nL\N^o{t) |ni) = nAle ^^5 
As\ {t) As) = nA\e~^^^ . 


iV(0)J \ns) , 

(49b) 

(50a) 

(50b) 


So, hnally we see that we can detect a fraction of order A\ of the other flavor in multiparticle 
short- or long-lived neutral kaon states. 


V. CONCLUSIONS 

We have shown that it is possible to formulate the consistent and probability-preserving 
description of the UP-symmetry-violating evolution of a system containing any number of 
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decaying particles. This has been done within the framework of quantum mechanics of open 
systems based on the approach developed in 1^, To achieve this aim we have considered 


master equations built up from creation and annihilation operators which generate dynamical 
semigroups that can describe the exponential decay and flavour oscillations for a system 
of many particles. It should be noted that this dynamical semigroup was used for the 
description of the entire system of unstable particles, and not only for decoherence effects, 
as it was done in, e.g., 00. 


We have used the fact that the decay of the particle can be regarded as a Markov process 
(the probability of decay of a particle is constant, so it does not depend on its history) 
and we have solved explicitly the Kossakowski-Lindblad master equation for a system of 
particles with violated CP symmetry. To avoid working with inhnite numbers of density 
matrix elements we performed calculations in the Heisenberg picture of quantum mechanics. 
This allowed us to deal with a low-dimensional system of linear differential equations for 
evolution of any observable bilinear in creation and annihilation operators. The choice of a 
concrete observable is then reduced to the proper choice of initial conditions for the system 
of differential equations. To show explicitly the effectiveness of the introduced approach we 
have calculated the evolution of the observables which are most interesting from the physical 
point of view. We have also found the evolution of mean values of these observables as well 
as their lowest order difference with the CP-preserved values. 

It seems to us that the presented analysis of time evolution of neutral kaons could be 
applied to the description of EPR states. Of course one must construct such a state by 
labeling annihilation and creation operators either by a discrete index, as it is done in 
quantum optics, or by a continuous parameter, as it is done in quantum held theory (beware 
that in the latter case annihilation and creation operators are actually operator valued 
distributions). 

Moreover, we think that our approach would be interesting and helpful also in analysis 
of mechanical systems and electric circuits. This follows from the fact that there is a strong 
analogy with nonreversible classical mechanics and quantum systems with CP violation j^ . 

Finally, we would like to point out that all the presented results are also valid for neutral 
B mesons after appropriate change of notation and values of physical quantities. This follows 
from the fact that B mesons evolve according to the same scheme as kaons. 
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